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Abstract 

The metric f(R) theories of gravity are generalized to five-dimensional spacetimes. 
By assuming a hypersurface-orthogonal Killing vector field representing the compact fifth 
dimension, the five-dimensional theories are reduced to their four-dimensional formalism. 
Then we study the cosmology of a special class of f(R) = aR m models in a spatially 
q-( flat FRW spacetime. It is shown that the parameter m can be constrained to a certain 

range by the current observed deceleration parameter, and its lower bound corresponds 
to the Kaluza-Klein theory. It turns out that both expansion and contraction of the 
extra dimension may prescribe the smooth transition from the deceleration era to the 
acceleration era in the recent past as well as an accelerated scenario for the present universe. 
Hence five-dimensional f(R) gravity can naturally account for the present accelerated 
expansion of the universe. Moreover, the models predict a transition from acceleration to 
deceleration in the future, followed by a cosmic recollapse within finite time. This differs 
• ■ from the prediction of the five-dimensional Brans-Dicke theory but is in consistent with a 



recent prediction based on loop quantum cosmology. 
PACS numbers: 04.50.-h, 98.80.Es 
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1 Introduction 

Since 1998, a series of independent observations, including type la supernova [US], weak lens jl], 
cosmic microwave background anisotropy [5], large scale structure [6j[7], baryon oscillation [EJIS], 
etc., have confirmed that our universe is undergoing a period of accelerated expansion. This 
result conflicts fiercely with the prediction of general relativity, and therefore triggers a "golden 
age" of cosmology since researchers in cosmology, astrophysics, general relativity, particle physics 
are all involved in seeking for viable models to explain this phenomena. Within the frame of 
general relativity, the cosmic speed-up can be viewed as an indication that the present universe 
is dominated by certain mysterious fluid with large negative pressure, called "dark energy". 
However, such simple explanations could hardly be satisfactory. The simplest ACDM model, 
which fits the observation data best so far, assumes a cosmological constant to be responsible 
for the dark content. But the observed value of the constant A is unnaturally smaller than any 
estimation by tens of orders [iTilfrT] . Other dynamical dark energy models then waive evolving 
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scalar fields to replace and resemble the cosmological constant (for review see p2]). These early 
attempts provide us with initial understanding of the complexity of the problem, though they 
chiefly serve for empirical fittings with comparatively poor theoretical motivation, and none of 
such models turn out to be problem- free. 

In view of the challenge and the fact that gravity is the only dominant long-range interaction 
insofar as we know, it is reasonable to consider the possibility that we have not fully understood 
the character of gravity on a cosmological scale. Therefore, the exploration of alternative theories 
of gravity are proposed. One of such examples is the Kaluza-Klein (KK) theory, which was 
initiated by the motivation of unifying the gravitation field and the electromagnetic field in a 
five- dimensional (5-D) metric. Since the fifth dimension is supposed to be a compact S 1 circle 
with an extremely tiny radius, it would actually yield no observable effect. Graceful though, the 
original version of KK theory fails in passing the Solar-System experiment [13] . But recently a 
new model of modified (non-compact) Kaluza-Klein cosmology was investigated in [13], where 
the universe in turns inflates, decelerates, and then accelerates in respectively early times, 
radiation dominated era and matter dominated era. This result meets the observational facts 
roughly. But there still exist problems such as that it leaves no decelerated period in the matter- 
dominant era for the formation of large scale structure. Besides, in [15] Brans-Dicke (BD) theory 
of gravity was introduced to match the Mach's principle. It is found in [16] that 5-D BD theory 
can naturally predict the cosmic acceleration without the requirement of a time-varying BD 
parameter u or a fabricated potential in its 4-D counterpart [TTtfTS]. Thus, it lends confidence 
and motivation to the effort of extending gravitational theories in 4-D to 5-D space-time. Also, 
note that in [16] the allowed range for parameters u — n is considerably wide, where u is the BD 
parameter and n represents the interaction between the extra dimension and other spatial ones. 
We may suspect that further information concerning the different choices of u> — n are masked in 
the formalism of 5-D BD theory. Given the argument for the equivalence (by evolutional effect) 
of the f(R) theory of gravity to certain special cases of BD theory up to a potential term [19], 
we are led to consider the 5-D f(R) theories of gravity for such information. It is worth noting 
that the scale of the compact extra dimension of 5-D gravity theories, which we are considering, 
is constrained to be less than 50 fim by current experiments [20H22]. 

The f(R) theories of gravity extends the Hilbert-Einstein action in general relativity to 



where f(R) is an arbitrary function of the curvature scalar R. The reason for such exploration in 
4-D spacetime is twofold. On the one hand, recent research has shown that a small correction 
to the Hilbert-Einstein action by adding an inverse term of R would lead to acceleration of 
the universe [23], and a large variety of models are proposed since then (for review paper, 
see [24J). However, few of them turn out to be without problem and most of them admit fine- 
tunings so as to be consistent with experiments. On the other hand, it is highly possible that 
future quantum theory of gravity, such as string theory, loop quantum gravity, etc., would bring 
about modifications to the action for classical general relativity, just like that general relativity 
has introduced corrections to the Newton's gravitation theory. Therefore, this effort would 
potentially match the low energy effective theory of quantum gravity and therefore accelerate 
its birth. 

With above motivations, we thought that modified gravity without fine-tuning should better 
be constructed in higher dimensions. It is shown in this paper that values of the parameters in 
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the 5-D theory are all with reasonable deductions and interpretations. In 5-D spacetimes with 
a hypersurface-orthogonal Killing vector field, the 5-D f(R) theories of gravity are reduced to 
the sensible 4-D world in section |2j Here we have adopted KK's idea that the fifth dimension 
is a small unobservable compact ring S 1 , so that a Killing vector field would arise naturally in 
the low energy environment. Further, in view of that in KK's theory one can find the sensible 

4- D universe to be expanding (without acceleration) as a result of the contraction of the extra 
dimension, we assume that there is certain interaction between the fifth dimension and other 
ones at present, and take their relationship to be the general power-law form. Then the reduced 
theories are carefully studied in the k = case of the FRW metric in section [3J We consider a 
special class of models with f(R) = aR m and present the numerical simulation of the evolution 
history of the matter dominated universe in section HJ demonstrating its consistency with the 
observation. Finally, we end with discussions in section [5] on the similarity and difference 
between the 5-D f(R) theories and the 5-D BD theory, as well as further indications from our 

5- D models. 



2 Killing reduction of the 5-D f(R) theory 

We start with the action in 5-D space-time: 

Here k = 8nG^/c 4 , where denotes the gravitational constant in the 5-D spacetime, and 
Sm represents the matter term in the total action S. Variation of (T5]) with respect to the metric 
g ab gives 

f\R)R ab - l -g ab f(R) - (V a V 6 - g ab V c V c )f\R) = nT ab) (3) 

where T ab = "^^fpF; an d f'(R) = df(R)/dR. Note that we have employed the abstract 
index [25] and R, R ab , T ab represent quantities in the 5-D spacetime. Contracting equation 
with g ab , we obtain the dynamical equation for the scalar field f'(R)- 

Waf'iR) = \[*r- Rf(R) + lf(R))- (4) 

Since here g ab g ab =5 instead of 4, the dynamical equation of f'(R) differs from its usual 4-D 
counterpart. 

Next we consider the structure of the 5-D spacetime and its relation to our sensible 4-D 
world. First, as mentioned in the introduction, we assume that the 5-D spacetime possesses a 
Killing vector field £ a which represents the fifth dimension and is everywhere space-like. Recall 
that the Killing reduction of 4-D spacetime is studied by Geroch [26] and is further extended 
to the 5-D spacetime by Yang et al. [27] . Therefore, following these works we introduce the 5-D 
metric as 

9ab = Kb + A^faft, (5) 

where h ab is the metric in the usual 4-D universe and A = £ a £ a . If £ a is not hypersurface- 
orthogonal, it would be related to electromagnetic 4-potential in the reduced 4-D theory, as 
demonstrated by KK theory. Since we are chiefly concerned with the cosmological effect of the 
reduced model, our discussion is restricted to the case where £ a is hypersurface-orthogonal for 
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convenience. Without losing generality, a coordinate system can be chosen as {x^x 5 },/! = 
0, 1, 2, 3, with {-^) a = i a - Then the line element of g ab reads ds 2 = g^dx^dx" + Xdx 5 dx 5 , while 
Ricci tensors and the dynamical equation of A in 5-D and 4-D spacetime bear the following 
relation [27] : 

,(4) 



R 



„, - ^X~ 1 D a D b X - l -\~\D a \)D h \ + h™h n b R mn , 

D a D a X = ~\-\D a \)D a \ - 2R ab Ce- 
Here D a denotes the covariant derivative on 4-D spacetime and is defined as 
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D a T, 
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satisfying all the conditions for a derivative operator. On the other hand, the stress-energy 
tensor in ([3]) is regarded as a perfect fluid in 5-D spacetime with the expression [16] 



r(5) 

" ab 



L- L \-2[(p + P)U a U b + Pg, 



ab 



Ti 4) + PXUb 



(8) 



where p and P are the 4-D energy density and the hydrostatic pressure respectively, and L 
is a constant representing the coordinate scale of the fifth dimension. This extension of the 
stress-energy tensor of the 4-D perfect fluid to a 5-D one can be understood in the following 
way: because the fifth dimension is compact and attached to every point of the 4-D space-time, 
p and P experienced by a 4-D observer should be an integrated effect throughout the compact 
ring. If we further expect that the fluid distributes homogeneously and does not travel along 
the fifth dimension, it is clear that p = J p (5) \ 1/2 dx 5 = p (5) A 1/2 L, P = P (5) A 1/2 L, and hence 
equation (jSJ) is obtained. Combining (j3J), (J1J with (JBJ, (jTJ), and (JSJ) , straightforward calculations 
lead to the 4-D field equation: 



G 



(4) 
ab 



8nG A-3 
c 4 f'(R) 



4- h 
L ab + h ab 



f(R) - Rf'(R) 

2f'(R) 



+ h-\D a D b - h ab D c D c )X - h- 2 [(D a X)D b X - h ab (D c X)D c X] 



f'(R) 



(D a D b - h ab D c D c )f'(R) 



2f'(R) 



X- 1 h ab (D c X)D c f(R), 



and the dynamical equations of A and f'(R)' 
8nG A5 



D a D a X 



c 4 f'(R) \2 2 P 
X 



l -X-\D a X)D a X - jl^(D a X)DJ'(R) 



-f(R) - -Rf'(R) 
f'(R) L4 2 J v ' 



D a D a f{R) = ^-^(T^+p) - l -X-\D*X)D a f{R) + \ 



-J(R) - Rf'(R) 



where G = G^L 1 represents the usual 4-D gravitational constant [28J 
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3 FRW cosmology of the reduced f(R) gravity 



In this section, we will study the cosmological predictions of the reduced f(R) gravity. The met- 
ric of a spatially isotropic and homogeneous 4-D spacetime is the Friedman-Robertson- Walker 
(FRW) metric with three possible structures of the space. As suggested by the observation [29J, 
we only handle the spatially flat case where the 4-D line element reads 

3 

ds 2 = -dt 2 + a 2 (t) dx 2 . (12) 
i=i 

Then the two components of the field equation (Q are: 

*( h \* - « n _ 3hX _ o* f"( R )n _ l* f"(R)R _ f(R)-Rf'(R) 
6 \a)- MG f>{R) P 2aX a f'(R) 2 X f'(R) 2f(R) {U) 

and 

2 a + a 2 \ = SttG Vj p | a A 1 f xY | lX f"(R)R | 2 <x f"{R)R 



a a?) c 2 f'(R) aX 4 I Ay 2 A f'(R) a f'(R) 

X f"'(R)R 2 + f"(R)R f{R)-Rf'(R) 

X f'(R) 2f'(R) ' 1 1 

The dynamical equations of A and f'(R) are respectively: 

X *„r x ~ 2 (f>\ ,",'AV \ f"(R)R f - 2Rf'(R) 

and 

f" {R)R 2 + f"(R)R 8nG ±i_ (P_P\_ 3 * f'( R ) R 1 * f"(R)R lf( R ) ~ R f'( R ) 



f'(R) f'(R) V4 c 2 J a f'(R) 2 X f'(R) Af'(R) 

(16) 

Here an over-dot denotes the derivative with respect to the proper time t of the isotropic ob- 
server, and we have assumed that the scalar field A depends only on t. Note that the conservation 
equation of the 5-D stress-energy tensor, V^T^q = 0, gives: 

O, i\ 1XP , , 

" + 3 a (p + ^ + 2^ =0 - (17) 

which can also be obtained directly from equations ( 1131) - ( TT6l) . For the common baryonic matter 
content of the present universe, the pressure is negligible compared to the energy density. Thus, 
we set P/c 2 = by approximation in the above equation and get p = p (a /a) 3 , where p is the 
current observed energy density of the luminary matter. 

By now, we have been dealing with the general form of f(R) gravity. Next we restrict our 
discussion to a specific class of f(R) models, f(R) = aR m , to further study the evolutional 
characteristics of the theory. It has been shown that this choice in 4-D spacetime, with m/1 
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and the common baryonic matter, could fit the Hubble diagram of Type la supernovae without 
need of dark energy |30j. We denote f'(R) = amR 771 ^ 1 = and obtain 



i i 

m-l rh I rh \ m-1 



R=( — \ , f(R) = - 4 

where a is a dimensional constant. With the above analysis, we find the three cosmological 
evolution equations for numerical simulations from the combination of equations (fl3l) - (fl6l) : 

2 ' ' l 

a ( a\ dA d0 1X6 ( 3 \\ ( 6 A™- 1 i ,3 /a \ 3 

a \a ) a X aq> 2 A \ 8m 4 / y am / 4 V a J 

A^A^a] -A?-U-2j(£J +8 ^ A ^ >U ' (20) 

7 = -3-7 " -Tli ~ i- - ~ )( + SnGX-h-^-po . (21) 

a0 2A0 \8m 4 J \am J Y A h \ a J K J 

4 Numerical simulations 

Let us now find the natural initial values of „o, do, Ao, Ao, 0o, 0o for the numerical simulation. 
Firstly, ao, Ao themselves have no direct physical meaning and therefore can be simply fixed as 
1 (with no dimension). Then we can directly determine the value of do by the present value 
of Hubble parameter H = (a/a) to . Secondly, we adopt the dynamical compactification idea of 
KK cosmology that the contraction of the extra dimension would result in the expansion of the 
remaining dimensions [31], and assume 

fl 3 A n/2 = constant ^ (22) 

where for a particular f(R) theory with given m and a in Eq. (1181) . the free parameter n can 
be constrained by the observation. Therefore we have (A/A)t = 6Hq/u. Finally, comparing the 
matter term in equation ([_]) to the field equation of general relativity, we expect X~ l l 2 S~ l ~ 1, 

at least for the present period; thus <f>o = 1 (with no dimension either) and further 
In summary, the initial conditions are: 



to 



1 _ 

2 I A 



to 

«o = A = 0o = 1 



do = Ho, Ao = —~H , 0o — -Hq. 



6 U X — 3 u 

(23) 



Then we estimate the value of a in equations (I19p - (12 ip and determine the allowed range for 
the pair of parameters m and n. A calculation of the 5-D curvature scalar shows 

/a d 2 \ A 1 ( X\ dA 

R = e {- a + 7<) + -x + 2{x) +3 ~x- (24) 

Inserting (ITS]) . fl20|) into the expression of R and noting that the deceleration parameter q = — ™ 
is an observable with its present value denoted as go, we have 



24m 



am I 1 + 2m 



<?o + -~r + 



9 8nGp 



n 2 12H% 



Hi (25) 
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Figure 1: Allowed range for the values of m and n with go £ ( — 0.47, —0.67) [32]. The horizontal 
axis stands for different value of n, while the vertical axis represents m. 



Substituting f )25|) into ffl~9]) . we obtain the relationship between parameters m and n as follows: 

5 - 4g + f 
2(1 - 2g ) + f + £gf 



where G = 6.67x Kr 11 ^- 1 ™ 3 ^- 2 , p = (3.8±0.2) x 10- 28 kg m - 3 , H = (2.3±0.1) x Kr 18 s" 
and therefore 87r Sf° w 0.12. One necessary and sufficient condition for the present accelerated 
expansion of the universe is that the deceleration parameter q < 0. Specifically, the allowed 
range for g according to the present observation is go = (—0.57 ±0.10) [32]. Here we apply this 
criteria to determine the allowed range for m and n, as depicted in Fig. [TJ From the figure, 
it is obvious that there do exist suitable parameters m and n, that is, suitable formalisms of 
f(R) = aR m and corresponding relationships between the extra dimension and the other ones, 
such that the present cosmic acceleration can be explained without the need of dark energy. 
Besides, Fig. 1 also reveals that both positive and negative values of n is allowed, indicating 
that the idea of dynamical compactification can be extended to that both the expansion and 
the contraction of the extra dimension can be responsible for the accelerated expansion of the 
other dimensions. Besides, we know from equation ( 125]) that 



1 

a = — 
m 



21m ( 9 8nG Po \ 



1 + 2m V n 2 \2H% J 



(27) 



Using Eqs. (1261) and (1271) . we depict the a — n relation in Fig. [2] with the range spanned 
by g , from which we find a increases rapidly as |n| increases. Further, one could find from 
Fig. [T]or equation (12"B"|) that the lower bound of m is 1, corresponding to a = 1, and therefore 
f(R) = R, which is exactly the KK case. One subtlety concerning this case is that Fig. [T]as well 
as equation (126]) shows that n = for m = 1, and therefore equation (122]) apparently indicates 
Ho = a = 0, contradicting with the present observations. Thus one might think that the KK 
theory is inconsistent with the observation. However, this conclusion is misleading, because 
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Figure 2: The relation between n and the possible range for a. The range is spanned by the 
uncertainty of the deceleration parameter q . 



from the inverse-solving procedure in equation (TT8|) . we have already confined our discussion to 
m^l cases. On the other hand, the asymptotic value of m, as n approaches infinity, is 

r 5-4g 

lim m = —, : , 28 

n^oo 2(1-2%) + 0.12' v ; 

from which we find the upper bound for m to be 1.72. Hence, the possible range for m is 
narrowed down as 1 < m < 1.72. 

With equations ( TT9l) - (12 ID . initial values ( 1231) and ( 1251) . the numerical simulation of the 
evolution of the scale factor a(t), the scalar field X(t),<f>(t) and the deceleration parameter q(t) 
can be conducted. As an example, we choose m = 1.5, n = 5.4. Then from equation (T26|) we 

get q = —0.6. The corresponding value of (— ) m_1 is approximately 17.3-f/g- It is useful to 
define the effective gravitational "constant" from Eqs. (JSj) - f TTTT) as 



G = G\' 1/2 f'(Ry 1 = G\- 1,2 <p-\ (29) 

Note that G is actually an evolving scalar. Then the evolution characters of a(t), q(t), X(t), 4>(t) 
and G are illustrated in Figs. |3]- [6] respectively. From these figures, it is clear that the present 
deceleration parameter q = —0.6 is in consistent with observation and the above analysis. 
Moreover, we note the following: 

1. q(t) rolls from a positive value to a negative one smoothly in the recent past. Specifically, if 
the age of the universe is 13.5Gyr [32], the universe turns from deceleration to acceleration 
at t — (13.5 — 1.2)Gyr. Thus, the "coincidence" problem in dark energy can be addressed 
in this case. 

2. After q(t) reaches the minimal value, it rolls back again and becomes positive at t — 
(13.5 + 13.4)Gyr. Thus, the universe will become decelerating in the future rather than 
be endlessly accelerating. 
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Figure 3: The evolution of a(t) (solid 
line) and q(t) (dashed line) with m = 
1.5, n = 5.4. Here t=0 corresponds to 
today, and the present value of the de- 
celeration parameter is go — —0.6. 
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Figure 4: The evolution of A(t) with 
m = 1.5, n = 5.4. Note that \(t) is 
decreasing at t = 0. 
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Figure 7: Numerical simulation for a longer period of cosmic 
evolution. The deceleration parameter q(t) becomes divergent 
at t — (13.5 + 47.6)Gyr(dashed line), as a(t) (solid line in the 
subfigure) approaches a constant value. 



3. Numerical simulation for a long period of cosmic evolution shows that the value of q(t) 
becomes divergent at t = (13.5 + 47.6)Gyr, as depicted in Fig. [71 Since a(t) > and its 
evolution is generally slow, the only reason for the divergence of q(t) is d — > 0. In other 
words, the universe would come to a static point in the finite future. Actually in this 
period, numerical simulation shows that all quantities, such as X(t) and <f)(t), still remain 
well-behaved. After this static point, a{t) starts to decrease for a certain period before 
numerical simulation fails. Thus, around the static point of a(t), numerical simulation is 
still effective. 

4. X(t) decreases from a large value in the past and increases slowly at present, while 4>(t) 
increases to a maximum value for the present period and decreases slowly in the future. 
After the static point of a(t), nonetheless, the decrease of <j)(t) becomes faster and therefore 
it would approach in the finite future when the numerical simulation fails. Actually, since 
<p appears in the denominators of Eqs. f|T9|) - f )2T|) . the analysis from this set of equations 
is no longer applicable for the evolution character beyond the = point. Instead, one 
should take advantage of original equation (EJ) (or Eqs. §S§ - ffTTj) by multiplying f'(R) 
on both sides). We leave the details about the static point and the evolution beyond to 
future investigations. 

5. The evolution of G is generally slow in the future. The present value of G is maximal 
with an almost zero evolving speed. Thus, it is also in consistent with the observation 
data that the present gravitational constant has a negligible evolving speed. 

On the other hand, as indicated by Fig. [TJ the expansion of the extra dimension would 
also bring about the acceleration of other dimensions. If we take m = 1.5, n = —5.4 instead of 
n = +5.4, the numerical simulation will give similar evolution features of a(t), q(t) and G(t) with 
the present value of deceleration parameter q = —0.6, which are consistent with observations. In 
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this case, the universe turned from deceleration to acceleration at t = (13.5 — 6.1)Gyr, which was 
4.9Gyr earlier than that in the n = 5.4 case. Correspondingly, the future transition point of the 
universe from acceleration to deceleration is also moved ahead to t — (13.5 + 8.0)Gyr, rendering 
that the total acceleration period of the universe changes only moderately from 14.6Gyr to 
14.1Gyr. Thus, the n < case is very similar to the n > case, except that the cosmic 
transition from deceleration to acceleration and the inverse occur both ahead of time. The 
reason for such temporal translation can be viewed in the following way. As shown before, 
when 1 < m < 1.72, the future asymptotic evolution of A should be increasing. Therefore, if the 
initial (present) evolution rate for A is set to be increasing at present rather than be decreasing, 
there should be a temporal translation with respect to the evolution of A and therefore to other 
quantities that can be affected. 

5 Discussions 

To clarify the argument for the equivalence of the f(R) theory of gravity to certain special 
cases of BD theory up to a potential term [19], we now compare the cosmic evolution of 5-D 
f(R) models developed in previous sections to the 5-D BD theory [IS] . From equation (jSJ), it 
is noticeable that the potential term is ^-^ttM — , and similar terms also exist in equations 

([TO]) and pip . These terms correspond to the (^;) m ~ 1 terms in f|T9|) - (T2~T]) . and would be 
shown to play significant roles in predicting the future evolution of the universe. There are 
three independent functions in the 5-D BD theory [16], namely a(t), X(t) and 4>(t), where 
4>(t) represents the core idea of BD theory on an evolving gravitation "constant" due to the 
interaction between the local gravitational fields with the faraway matter in the universe. We 
note that q(t) in this case (see Fig. 2 in [16]) evolves to a constant in the future. This means that 
the universe would accelerate permanently in the future until its energy density decreases to 
be extremely low so that certain quantum effects may become significant, such as the situation 
in [33]. It can also be read from Figs. 3 - 5 in [16] that the asymptotic values of ^ and ^ are 
both zero, as further illustrated in Fig. |HJ 

In contrast, in our 5-D f(R) theories, Fig. [5] has demonstrated that, while <ft evolves to a 

large value presently and remains positive in the future, the (— ) m ~ 1 terms render the evolution 
of A(t), 4>{t) to be more complicated. Specifically, the corresponding evolution characteristics of 
^, j, 4 are illustrated in Fig. 0with specific parameters m = 1.5, n = 5.4. As G = GA 1 / 2 ^ 1 is 

an observable, it is interesting to show how (— ) m_1 terms affect its evolution. Since m > 1, 

' ° V am I 

l 

go < 0, from equation fl25l) we have (^;) m_1 > 0. Moreover, we do not consider any solution 

with a negative because of its anti-gravity character. Thus, the coefficient of the (^) m_1 
term would determine the evolutional inclination of A and 0. Specifically, because m = 0.5 and 
m = 2.5 correspond to — |) = and (^ — |J = respectively, the qualitative features of 
the future evolution of A and 6 follows: 

• If m < 0.5, A 4, (f) I 

• If m > 2.5, A|, 0t 

• If 0.5 < m < 2.5, At, <H 
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Figure 8: The evolution rate of ~ (solid Figure 9: The evolution of - (solid 
line), j (dashed line) and | (dash- line), ~ (dashed line) and | (dash- 
dotted line) in the 5-D BD theory, with dotted line) in 5-D f(R) gravity with 
parameters u = 1.2, n = 3 in [16]. parameters m = 1.5, n = 5.4. 



In view of the constraint illustrated by Fig. [H only the last situation is allowed, where the 
evolution of A and <fi counterbalances each other and therefore the evolution of G is generally 
slow. 

In summary, we present the Killing reduction of 5-D f(R) theories of gravity to the 4-D 
sensible world. Then we study its cosmological implication by assuming the spatial homogeneity 
and isotropy, namely, the Friedman- Robert son- Walker metric. With the illustration of a specific 
example it is found that the theory is consistent with the present observations in a variety of 
aspects, including the recent cross from cosmic deceleration to acceleration, its prescription of 
present cosmic speed-up, the negligible evolving speed of the effective gravitational constant, and 
etc.. In contrast to the KK cosmology, it is worth mentioning that in 5-D f(R) theories of gravity, 
both expansion and contraction of the extra dimension could result in the present accelerated 
expansion of other spatial dimensions. Even for a simple and generic class of f(R) = aR m 
models, the 5-D f(R) theories of gravity do not need unreasonable or ill- initiated fine-tuning 
of certain parameters. Hence it is reasonable to infer that the present accelerated expansion of 
spatial dimensions could be certain basic character of the 5-D spacetimes. 

Finally, it is pointed out in Ref. |33j that in the large scale limit, quantum gravity effect would 
come to be vital and bring about higher-order quantum corrections to the Friedman equation to 
prescribe a static point in the finite future, followed by the recollapse of the universe. Note that 
the static point in the finite future is also a significant character of the 5-D f(R) gravity. Thus, it 
is of interest to explore whether there are some relation between the 5-D theories of f(R) = aR m 
models and the quantum gravity theory applied in [33]. In addition, the future recollapse also 
differs from the prediction of endless acceleration scenario in the 5-D BD cosmology [16J. Hence, 
5-D f(R) gravity is not simply equivalent to 5-D BD theory, since the effect of potential terms 
are shown to be vital enough to prescribe different cosmological scenarios in the distant future. 
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